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ABSTRACT 

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations 
were carried out using the HEterodyne Receiver Array (HERA) at the 30m IRAM telescope in the CO(2-l ) line achieving a resolution 
of 12" x2.6 km s -1 , enabling individual Giant Molecular Clouds (GMCs) to be resolved. The observed region is 650 square arcminutes 
mainly along the major axis and out to a radius of 8.5 kpc, and covers entirely the 2' x 40' radial strip observed with the HIFI and 
PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK 
at 2.6 km s -1 velocity resolution. The CO(2-l) luminosity of the observed region is 1.7 ±0.1 x 10 7 K km s~'pc 2 and is estimated to 
be 2.8 ± 0.3 x 10 7 K km s~'pc 2 for the entire galaxy, corresponding to PL masses of 1.9 x 10 s Mq and 3.3 x 10 s Mq respectively 
(including He), calculated with a A , (H 2 )/Ico<i-0) twice the Galactic value due to the half-solar metallicity of M33. Hi 21 cm VLA 
archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in the CAS A software package, 
yielding a series of datacubes with resolutions ranging from 5" to 25". The Hi mass within a radius of 8.5 kpc is estimated to be 
1.4 x 10 9 M©. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9 ± 0.1 kpc whereas 
the atomic gas surface density is constant at S H i = 6 ± 2 Mqpc~ 2 deprojected to face-on. For a A , (H 2 )/Ico(i-0) conversion factor twice 
that of the Milky Way, the central kiloparsec Hi surface density is Zh, = 8.5 ± 0.2 Mqpc~ 2 . The star formation rate per unit molecular 
gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to H„ and FIR brightness, is 
constant with radius. The SFE, with a A'(H2)/Ico<i-0) factor twice galactic, appears 2-4 times greater than of large spiral galaxies. A 
morphological comparison of molecular and atomic gas with tracers of star formation is presented showing good agreement between 
these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated 
more than 8 kpc from the galaxy center, are presented. 

Key words. Galaxies: Individual: M33 - Galaxies: Local Group - Galaxies: evolution - Galaxies: ISM - ISM: Clouds - Stars: 
Formation 



1. Introduction 

The Local Group galaxies span a broad range in mass, luminos- 
ity, morphology, and metallicity. Two large spirals (the Milky 
Way and M31) are the centers of two galaxy sub-groupings, 
each being surrounded by a large number of dwarf galaxies. In 
addition, M3 1 — the Andromeda Galaxy — has a small spiral 
companion, M33 (the Triangulum Galaxy); their separation is 
approximately 15 degrees, corre sponding to 200 kpc (assuming 
a common distance of 840 kpc; iGalleti et alj d2004l) ). Gaseous 
streams are observed between them, indicating tidal interac- 



tion dPutman et al.ll2009h . M33 provides a means of observing 
a galaxy morphologically similar to our own but with a mass 
only a tenth of the Milky Way and factor two low er metallicity 
(iRosolowskv & Simonll2008l: iMagrini et alJl2009l) . Further evi- 
dence for the difference between M33 and the Milky Way is the 
large gas fraction and blue stellar colors of the former relative to 
the latter. M33 thus represents an environment in which to study 
the interstellar medium (ISM) and star formation (SF) that can- 
not be replaced by Galactic observations and where individual 
GMCs can be resolved to probe their star formation activity. It 
may also be possible to apply what we learn by studying M33 to 
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Table 1. Adopted parameters for M33 



a(J2000) 
(5(J2000) 
Distance 

Optical Radius R25 

Inclination 

Position angle 

Central Oxygen abundance 



l h 33 m 50!9 
+30°39'39" 
840 kpc a 
30.8' b 
56° b 
22.5° b 

12 + log(0/H) = 8.4 c 



" iGalleti et alj j2004b 

* HYPERLEDA dPaturel et all2003h 

c iMagrini et all d2009h 



the physics of early-universe objects, which share many of the 
characteristics of M33. 

In this article we present sensitive and high-resolution map- 
ping observations of the CO / = 2 — » 1 transition in M33 in 
order to study the morphology and dynamics of the molecular 
component. The total mapped area covers 650 square arcmin- 
utes, mainly along the major axis of the galaxy. 

A 2' x 40' wide strip along the major axis (see Fig. [TJ 
was observed to a particularly low noise level of 25 mK at 
2.6 km s _1 velocity resolution to compare with the sensitive 
Cn Herschel/HIFI and Herschel/PACS spectroscopy observa- 
tions which wi ll be obtained as par t of the HERM33ES Herschel 
Key Program (iKramer et al 1 l2010h . While the most sensitive 
and among the highest r esolution, these a r e not the first maps 
of M33 in the CO lines. lEngargiola et al.l (120031) observed the 
whole of the inner disk (up to about 5k pc along the major a xis) 
with the BIMA array at 13" resolution: iHever et all d2004l) ob- 
served the inner disk and a small major axis s trip at 50" res- 
olution with FCPvAO; iRosolowskv et al.l (120071) combined the 
BIMA+FCRAO+NRO data to i mprove the sensitivity and res- 
olution of the previous maps; and lGardan et al.1 J2007) observed 
a rectangle at high sensitivity and 15" resolution extending 
from NGC 604 to the R25 radius. Table [2] summarizes the 
characteristics of previous mol ecular and atomic g as surveys 
in M33. This work extends the iGardan et al.1 (120071) work fur- 
ther North and to the South at higher resolutio n and sensi- 
tivity. Earlie r studies of GMC s in M 33 include IWilson et al.l 
dl997l) and IRosolowskv et al.l d2003h and studies similar to 
our own of other Local Group galaxies have been made b y 
e.g. iFukui et ail (120081) : llsrael et al.l (120031) : iLerov et all (120061) : 

lNietenetal.l(l2006l) . 

As men t ioned by iBlitz & Rosolowskvl (120061) for M33 and 



Lerov et al.l (120061) f or IC 10 and discussed more extensively by 



Gardan et all d2007l) . the Star Formation Rate (SFR) per unit H 2 
mass or Star Formation Efficiency (SFE = SFR/M(H2) in yr~') 
was found to be up to an order of magnitude higher in these small 
galaxies than in large spirals. This appears to be the case in dis- 
tant galaxies as well, g iven the factor 10 increase in (commov - 
ing) SFR density (e.g. iMadau et~allll996t Iwilkins et al.ll2008l) . 
Are there local universe analogs of these distant objects? Is M33 
one of them? 

One of the obvious questions is whether the H2 mass has 
not been underestimated in these subsolar metallicity objects. 
The articles using the data presented here and as part of the 
HERM33ES project will attempt to answer that issue clearly. 

Metallicities were lower in the past and H2 production is be- 
lieved to take place on grain surfaces - therefore H 1 to H2 con- 
version is expected to be less efficient in low metallicity systems 




Fig. 1. The Local Group galaxy M33. This color image shows 
the GALEX FUV data which trace young stars and dust in the 
disk through attenuation. Overlaid are the observed fields (a) us- 
ing the IRAM 30m Pico Veleta telescope in the CO(2-l) line 
(thick white outline) and (b) by HIFI/PACS instruments as part 
of the HERM33ES Herschel Key Program (thin yellow stripe). 



dKrumholz et al.l2009|). and this is indeed observed in many sys- 
tems (e.g. lLerov et al.ll2007l) . The conversion of H 1 to FL in the 
intermediate redshift systems would have to be much more effi- 
cient than today, generating unusually high molecular fractions 
in these distant galaxies (contrary to expectations) with less stel- 
lar gravity, in order to have a similar efficiency in converting H2 
to stars. 

Here, we will present new maps of CO and HI. In this first 
paper, we restrict ourselves to a study of the radial distribution of 
the molecular gas and infrared surface brightness, the molecular 
and atomic gas surface densities and the star formation efficiency 
and to a qualitative comparison between the maps of star forma- 
tion rate tracers, i.e. the dust maps, and the gas maps. In addition, 
we will discuss CO spectra in a few selected regions. A series of 
articles will follow, focussing on at least (/) cloud populations, 
life cycle, and mass spectrum (if) dynamics of the molecular gas 
and the role of spiral arms (Hi) diffuse CO emission, after sub- 
traction of the clouds identified and (iv) a more detailed com- 
parison of the star formation rate - gas surface density relation. 
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CO Survey 


Telescope 


Line Anj 


jular Spectral 


Area 




Sensitivity 








resolution resolution 






per channel 










(") ( km s ') (arcmin 2 ) 


(Kpc) 


(mR) 


This paper 


IRAM 30m 


CO(2- 


-1) 


1 _ 


Z.O 


643 


8 


15-30 


Gardan et al. (2007) 


IRAM 30m 


CO(2- 


-1) 


1 J 


Z.O 


253 


7 


15-30 


Rosolowskv et al. (2007) BIMA+NRO+FCRAO 


CO(l- 


-0) 




Z.O 


170 


2 


60 


Rosolowskv et al. (2007) 


BIMA+FCRAO 


CO(l- 


-0) 


13 


2 


900 


5.5 


240 


Hever et al. (2004) 


FCRAO 


CO(l- 


-0) 


45 


1 


900 


5.5 


53 


Engargiola et al. (2003) 


BIMA 


CO(l- 


-0) 


13 


2 


900 


5.5 


240 


HI Survey 


Telescope 


Angular 


Spectral 




Sensitivity 


Sensitivity 








resolution 


resolution 




per channel 


per channel 








(") 


(kms- 1 ) 


(kpc) 


(mJy/beam) 


(K) 


This paper 


VLA BCD 


5-25 a 


1.3 


9 


1.1- 


-2.8 a 


24-2.75 a 


Putman et al. (2009) 


Arecibo 


204 


5.15 


22 




55 


0.3 


Corbelli & Schneider (1997) 


Arecibo 


270 b 


4.1 


20 




c 


c 


Deul & van der Hulst (1987) 


WRST+Effelsberg 


100m 


12x24 


8.2 


9 






1.2 



" See Table[3]for details. 

* Undersampled hexagonal grid, the value corresponds to the grid spacing. 

The authors give an integrated intensity sensitivity value of 1-2 Jy km s or 5.5-1 1 K km s 



2. Observations and data reduction 

2. 1 . IRAM CO(2-1 ) observations 

The observa t ions presented here are a follow-up to the 
iGardan et"ai1 d2007l) mapping of a large part of M33. All data 
were taken with the 30 meter telescope run by the Institut de 
RadioAstronomie Millimetrique^ (IRAM) on Pico Veleta near 
Granada, Spain. The observations were carried o ut starting in 
Augus t 2008 and also include the data taken by IGardan et alJ 
(120071) . 

The IRAM HERA instrument (HEterodyne arRAy), com- 
posed of 9 dual -polarization 1.3 mm receivers dSchuster et al.l 
l200l . was used in On-The-Fly scanning mode. All data are 
presented in the main beam temperature scale and we have as- 
sumed forward and main beam efficiencies of T]f or = 0.90 and 
T}mb = 0.49 + 0.04 for the HERA ob servations, the sensitivity is 
then 9.6Jy/K dSchuster et al.ll2004l) . The scanning speed was 3 
arcseconds per second of time with dumps every second. Both 
the VESPA backend, at 1.25MHz resolution, and the WILMA 
correlator with 2MHz channel spacing, were used simultane- 
ously. Reference positions were observed roughly every three 
minutes and these were taken to be holes in the H i column den- 
sity outside of the M33 disk with no visible FIR emission. The 
average beam size over the nine pixels is 11.7". The array was 
tilted by 18.5 degrees, obtaining 7.6 arcsecond spacing between 
the individual pixel tracks, and then shifted in position by 3.8 
arcsec in order to obtain full sampling. This was first done for 
scanning parallel to RA, then parallel to Dec, the array was then 
rotated by 90 degrees and the procedure was repeated. The array 
rotation is to have different pixels cover the same regions to im- 
prove data homogeneity and reduce striping. For the HIFI low 
noise major axis strip, additional scans along and perpendicu- 
lar to the strip direction where acquired. The repetition of inde- 
pendent scans over the same areas also randomizes any pointing 



•' IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



errors. Thus, if there are pointing uncertainties of 2", these con- 
tribute negligibly to the overall beam size (through a negligible 
smearing). The part of M33 presented here comes from the ob- 
servation of about 5 million spectra. 

Reduction was carried out within GildasOCLASS and GREG 
software. The intensity of the CO lines is usually small so that 
emission is not seen in individual spectra, the baseline removal 
was therefore implemented using a windowing based on the H i 
data at 17" presented in section [2721 In a first step a linear base- 
line is fitted and subtracted from each individual spectrum not 
taking into account channels inside a window computed by find- 
ing the first channels at a 3<x level framing the peak channel of 
H i emission and adding a 40 km s _1 buffer on each side of this. 
A filtering step is then applied to remove spectra with anomalous 
noise. The actual noise is computed for each spectrum and com- 
pared to the theoretical noise computed from the system temper- 
ature, integration time and spectral resolution. The spectra with 
excess noise attributable to poor baselines are removed; this cor- 
responds to about 10% of the spectra taken. The remaining spec- 
tra are then convolved by a gaussian to obtain a regularly gridded 
cube of the desired resolution. We then fit and remove a third or- 
der polynomial to the spectra corresponding to each of the spatial 
pixels in the gridded cube this time reducing the velocity buffer 
described previously to 30kms~'. With this method we have 
made CO cubes with spacial resolutions of 12, 15 and 25" and 
with velocity resolutions of 2.6 km s -1 . 

We compute t he CO(2-l) integra ted intensity map using a 
masking method dGratier et al.1 1201 51) taking into account Hi 
data, we developed in order to filter out some of the noise 
present in the observations and increase the sensitivity to low 
intensity possibly diffuse CO emission. Previous masking meth- 
ods used masks created from spatially smoothed versions of the 
original CO data cubes to filter out re gions dominated by noise 
(lAdleret all [1991 iDigel et al.l [19961 lLoinard etal] fl999h . We 
use the 21 cm atomic hydrogen data at 17" x 17" x 1.27 km s _1 



b http://www.iram.fr/IRAMFR/GILDAS 
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resolution presented in section [2721 to achieve the same goal, the 
underlying hypothesis being that molecular gas is unlikely to be 
present for low enough values of N{H i) so the corresponding ve- 
locity channels can be discarded when computing the integrated 
intensity CO map. For each pixel of the H i cube, we create a 
binary mask keeping only the velocity range for each pixel cor- 
responding to a H i signal value above a defined threshold in the 
Hi signal. The integrated moment map for the CO(2-l) data is 
then computed summing only velocity channels included in the 
H i mask. Figure [2] shows the integrated intensity map for the 
12" cube using this method and Fig.|3]corresponding noise map 
at 12" computed from 30 channels free of emission from M33 
(from -80 to km s Vlsr). The result is an increased S/N 
ratio as the channels contributing only noise to the sum are no 
longer taken into account. The value of the H i signal threshold 
was chosen at 10 K which corresponds to 1.9<x where <x is the 
rms noise level in the 17" H i cube. We tested masking values be- 
tween 5 and 25 K (0.9 to 4.6<r) and the total CO intensity varied 
by only a few percent. Significantly above or below these values, 
CO signal is lost or more noise is included. At the 10 K bright- 
ness level, the Hi lines are broad and sample several channels 
to either side of the rotation curve. Comparing with other tech- 
niques of looking for CO, we found no evidence of missing CO 
emission by this technique. As an examp le, the "Lonely Clou d" 
(a = l h 34 m 16 s .7, 8 = +30°59'3", J2000) dGardan et alj|2007l) in 
an interarm region not near an H i maximum is included by this 
method. 

2.2. VLA H i mosaicing 

The H i observations are from archival VLA B, C, and D array 
data taken as part of projects AT206 and AT268 in 1997, 1998 
and 2001 . At the frequency of the H i line the VLA primary beam 
is 32'. The mosaic comprises 20 D array pointings on a nearly 
square grid separated by 16'and 6 B and 6 C pointings in a dia- 
mond shape covering the entire star forming disk. The primary 
beam centers of the B/C and D array mosaic are shown in Fig. [4] 
overlaid on the Spitzer 70/im map - the D array data contains 
more pointings extending further from the center of M33. Since 
our goal is compare the distribution and kinematics of the atomic 
and molecular gas at high resolution, we focus here on the stel- 
lar disk where VLA B-, C- and D-array data are available. Data 
calibration and imaging were carried out using the 2.4 version of 
the CAS A0 software package. The quasar 3C48 {a = l h 37 m 41 s :3, 
8 - +33°09'35") was used as flux, bandpass and phase calibra- 
tor (with an adopted flux of 16.04 Jy). The calibration of the vis- 
ibilities was done through the following tasks. First the correct 
flux was applied to the calibrator visibility model using setjy, 
then frequency bandpass and gain calibrations were computed 
using the bandpass and gaincal tasks. 

The calibration solutions were finally applied to the entire 
dataset with the use of the applycal task, continuum emission 
averaged on channels not including line emission from M33 was 
subtracted in the wv-space using the uvcontsub task. Flagging 
of bad data was done by hand. The imaging was achie ved using 
the mu lti-scale clean MSCLEAN algorithm described in lCornwelU 
J2008h . This method has several advantages over the classical 
clean algorithms that use point like sources to model the ob- 
served emission. By cleaning several scales simultaneously, it 
is able to efficiently recove r extended emission and elim inates 
problems of flux correction dJorsater & van Moorselll995l) . neg- 
ative bowls surrounding regions of extended emission, and the 

c http://casa.nrao.edu/ 
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Fig. 2. IRAM CO(2-l) integrated intensity map of the galaxy 
M33 at a resolution of 12". This map was made using the VLA 
Hi cube (17" resolution) as a mask, employing a 10 K cutoff 
per channel (see Sect 12.21 ). As a consequence, areas without HI 
emission above that level appear white. Grey ellipses are drawn 
every kpc using orientation parameters listed in Table Q] 



pedestal of low level uncleaned flux in the final residual map. 
An extensive comparison of multisc ale clean and its po int source 
counterpart has been carried out in iRich et al.l (|2008). A series 
of four cubes with angular resolutions of 5, 12, 17, and 25 arc- 
seconds and 1.27 km s channel width were computed. The 5" 
cube corresponds to the highest resolution achievable with the 
B array. The other cubes were obtained using gaussian tapers 
in the wv-plane of respectively 4900 and 3150 and 1900 meters 
FWHM. In order to min imize sidelobe s over the entire map, a 
robust weighting scheme iBriggsl (1995) was applied to the data 
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Fig. 3. CO noise map of the 12" cube computed from 30 chan- 
nels free of M33 emission, from -80 to km s _1 Vlsr- 
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Fig. 4. H i primary beam coverage of M33 overlaid on the 70^m 
Spitzer map. VLA D pointings are in red (solid lines) and B/C 
ones in yellow (dashed lines). 



Table 3. Beam and rms noise properties of the Hi 21 cm dat- 
acubes. 



Beam 


PA 




0~t 


(" x ") 


O 


(mJy/beam) 


(K) 


5.5x5.2 


-95.1 


1.1 


24 


12.0 x 11.6 


-31.8 


2.0 


9 


17.2 x 17.1 


-45.8 


2.5 


5.4 


25.9 x 24.2 


-74.8 


2.8 


2.75 



The rms noise, in mJy/beam and brightness temperature, was 
calculated over an ellipse of galactocentric radius 8.5 kpc after 
primary beam correction. 



using a robustness factor of 0.5. The rms noise for the cubes was 
computed over the channels without emission from M33 within 
a galactocentric radius of 8.5 kpc. Table|3]summarizes the beam 
sizes and noise properties of the different cubes. A zeroth mo- 
ment map (Fig. [5] with the CO 1 K kms~' intensity contour 
superposed) was computed at 17" resolution by masking out re- 
gions of the 17" cube where emission from the 25" cube was 
below 2.75K (lcr). The data presented here will be made avail- 
able throught the Centre de Donnees de Strasbourg 

(CDSfl 



d |http : // cdsweb . u- strasbg . fr/| 



3. Results 

3. 1 . Gas masses 

In this section we present measurements of total integrated flux 
and derived total molecular and atomic gas masses. 

3.1.1. Molecular Gas 

The CO(2-l) luminosity measured over our observed area 
is 1.7 + 0.1 x 10 7 Kkms _1 pc 2 . Our CO mapped region 
covers 28% of the total galaxy area up to 8.5 kpc but 
61 + 2% of the infrared emission (see details in Sect. 13.3b , 
The estimated CO(2-l) luminosity for the whole R25 disk 
is thus 2.8 + 0.3 x 10 7 Kkms^pc 2 . The molecular gas 



6 P. Gratier et al.: Molecular and Atomic Gas in the Local Group Galaxy M33 

Table 4. Comparison of results with previous studies of M33. 



CO surveys 


Molecular gas mass b 


SFRS 


Depletion time k 


Depletion time 1 




(M ) 


(M yr- 1 ) 


(10 8 yr) 


(10 8 yr) 


This study 


1.9xl0 8 (3.3xl0 8 ) c 


0.27 h 


7 


7 


Gardan et al. (2007) 


3.2 x 10 7d 


0.05-0.13' 


2-6 


1-3 


Rosolowskv et al. (2007) 


J.O X 1U 


U.4j J 


Q 
O 




Hever et al. (2004) 


3.5 x 10 


A A CI 

0.45 J 


O 

8 


1 O 

1-3 


bngargiolaet al. (2003) 


9 x 10 


0.24 


4 


1 


Corbelli (2003) a 


2.9 x 10 8f 


0.45 j 


6 




HI surveys 


Atomic gas mass 111 








This study 


1.4 x 10 9 








Putman et al. (2009) 


1.5 x 10 9n 








Corbelli & Schneider (19971 


2.6 x 10 9 ° 








Deul & van der Hulst (1987) 


1.1 x 10 9 









" FCRAO CO(1-0) data from lHever etai] J20041) . 

* All values converted to our adopted Xco = 4 x 10 20 cirT 2 /(K km s~') (twice Galactic). 

c The value in parentheses is the extrapolation to the whole galaxy for comparison with other surveys. See Sect. l3.1.T1 for details. 

d Original value used X co = 2 x 10 20 cirr 2 /(K km s" 1 ) (Galactic). 

e Original value used Xco = 3 x 10 20 cirT 2 /(K km s -1 ). 

- f Original value used X co = 2.8 x 10 20 cirr 2 /(K km s" 1 ). 

s Over the areas mapped in CO for each survey. 

h From lVerlev et al] d2009l) 

' From H„ with a|Kennicutt ( 1998]) calibration. 

i We give the value from[Engargiola et al. (2003) as the CO mapped areas are identical. 

* Value computed by dividing the first by the second column of this table. 
' Original value given by the cited papers. 

"' All values converted to our adopted distance D = 840 kpc. 

" Original value used a distance of 730 kpc. Inside a 3.2 x 10 20 cirT 2 contour. 

° Original value used a distance of 690 kpc. Over their whole mapped area extending beyond 8.5 kpc. 



mass is computed with the following hypotheses. Given that 
metallicity is half solar and the gradient in M33 is weak 
(12 + loz(OIH) = 8 .4 - 0.03jy, iRosolowskv & Simonll2008t 
Mag rini et al] 20091). we use a constant "standard" Milky Way 
factor ( Dickman et al.l Il986l) multiplied by a factor of two 



7Y(H 2 )/Ico(i-0) = 4x10 cmT 1 /{YL km s~'), implicitly assuming 
an inverse l inear depende nce bet ween the con version factor and 
metallicity (IWilsonl 19951) but seelisrael(l2000l) . Using the CO(l- 
0) observations from lRosolowskv et al] (120071) and our CO(2-l) 
observations we compute a CO(2-l)/C0(l-0) ratio of 0.8 in the 
central kiloparsec of M33. Due to the decrease of the excitation 
temperature with radius we can ex pect lower values of the line 
ratio for the ou ter parts as found by ISawada et al.l (12001 1) for the 
Milky Way and lBraine etail (1 19971) for NGC4414. With a linear 
gradient of the line ratio from 0.8 at the center of the galaxy to 
0.5 at 8.5 kpc, and taking into account the measured exponen- 
tial decrease of the CO emission with radius (see Sect 13.3b we 
find an average CO(2-1)/(1-0) line ratio of 0.73 for the whole 
galaxy. The computed mass includes helium by multiplying by a 
factor 1.37. The molecular gas mass for the whole galaxy is thus 
3.3 x 10 8 M Q . 



3.1.2. Atomic Gas 

The integrated H i flux density over a region corresponding to a 
galactocentric radius smaller than 8.5 kpc, slightly over the R25 
radius for M33, is 8.4 x 10 3 Jy kms" 1 . The atomic gas mass 



was derived from the H 1 integrate d intensity using a conve rsion 
factor 1.8 x 10 18 crrT 2 /(K km s 1 ) dRohlfs & Wilsonll 19961). The 
total H 1 mass over the same region is 1.4xl0 9 M^. lPutman et al.l 
(120091) using the Arecibo telescope find give 1 .5 x 10 9 Mq inside 
a very similar contour which they define as the star forming disk. 



3.2. Kinematics 

We obtain velocity fields for both datasets by computing the 
first moment of the emission. In the case of Hi, the 17" data 
was masked by the 25" cube with a threshold of 10K. For the 
CO data, we computed the first moment of the original non H 1 
masked datacube, masking out emission below 4<x. Both meth- 
ods are thus completely independent of each other. Fig.|6]shows 
the CO and Hi velocity fields at respectively 12" and 17". The 
excellent agreement between the two velocity fields is a further 
argument in favor of using the H 1 masking to obtain CO inte- 
grated intensity maps. 



3.3. Radial distribution 

The radial distribution of properties were derived taking az- 
imuthal averages over elliptical rings with inclinations of i=56° 
and position angle PA=22.5° (measured toward East from 
North). The quantities are therefore deprojected from the sky 
plane onto the plane of M33. The radial step is 500pc which 
corresponds to about 2'. The bottom part of Fig. [7] shows the 
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Fig. 6. CO (/e/f) and Hi (right) velocity fields at respectively 12" and 17". The CO velocity field was computed from the original 
non Hi masked cube, both velocity fields are thus independent. The contours on the right correspond to a constant S/N ratio of 10 
in the 12" CO integrated intensity. 



cumulative fractions of the area (solid line) and 70 yum surface 
brightness (dashed line) in our CO map (black area in the in- 
set) compared to the whole elliptical region of radius less than 
R (sum of grey and black areas of the inset). For galactocentric 
distances below 2 kpc our map samples more than 80% of the 
accessible disk at these radii both in area and 70 yum flux. The 
whole CO mapped area corresponds to 28% of the area up to 
8.5 kpc but more than 60% of the corresponding 70 /im flux, this 



is explained by the rapid decrease of infrared emission with ra- 
dius. The corresponding values for 24 and 160 /urn are similar to 
a few percent at each radius. Tab. [5] summarizes the source and 
resolution of the ancillary data used. 
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Fig. 5. VLA HI column density map (in color) of the galaxy 
M33 at a resolution of 17" overlaid with the IRAM CO(2-l) 
intensity contour at 1 K km s _1 . The region observed in CO(2-l) 
are marked with thin white lines. 



3.3.1. Surface brightness 

The top part of Fig. [7J presents the radial variations in average 
surface brightness of the CO and infrared data from Spitzer at 
24 yum, 70 jum and 160 /im, the solid lines correspond to av- 
erages taking into account only our CO mapped area and the 
dashed lines to the whole elliptical annuli for a given radius. 
For each tracer, the similarity between these two averages is an 
indication that the area mapped in CO is representative of the 
M33 disk in its entirety. The 4 datasets have similar exponen- 
tial decreases with radius i.e. F v (r) <x e~'/ L with a correspond- 
ing scale length L. Table [6] summarizes the scale lengths com- 
puted by least square fitting for CO and IR data, and Table [7J 
scalelength values from surveys presented in Table [2] The in- 
frared data scale length s are in agreement with the values in 
iTabatabaei et alJ d2007ah and lVerlev et al.l (120091) albeit slightly 
smaller as the IR tracers seem to drop marginally more rapidly 
in our CO mapped region (solid lines in Fig|7J compared to the 
whole galaxy (dashed lines in the same figure). The scale length 
of CO is higher than the ones of 24 fim and 70 pm and similar to 
the 160 /im. 
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Fig. 7. (Color online) Top: Radial distribution of the CO and IR 
surface brightness in units of MJy sr~' for infrared data and 
Kkms 1 for CO. For the dust tracers, the surface brightness 
is computed both on the region covered by our CO observa- 
tions (solid line) and on the whole elliptic annuli for each radius 
(dashed lines), showing that our CO sampled area is representa- 
tive of the galaxy for each radius. Bottom: Cumulative fractions 
of the area (solid line) and 70 pm. surface brightness (dashed 
line) in the region covered by our CO observations (black area 
in the lower left diagram) compared to the whole elliptic region 
of radius less than R (sum of grey and black areas at a given R). 



3.3.2. Mass surface densities 

We compute mass surface densities for atomic and molecular gas 
and study their radial distribution. For the molecular gas mass, 
we use the same hypotheses as in Sect l3. 1 . ll including a constant 
CO(2-1)/CO(1-0) factor. Taking a linearly decreasing value of 
this factor would modify only slightly the radial distribution, 
raising the most extreme outer points by a factor 1 .5 and low- 
ering the inner ones by a factor 1.1, at the price of added uncer- 
tainties that cannot be simply estimated. The top panel of Fig. [8] 
shows the mass surface density of atomic and molecular gas as a 
function of radius. The atomic gas surface density is constant at 
an average value of 6 + 2 M0pc~ 2 , the molecular gas mass sur- 
face density follows the same trend as the CO surface brightness 
presented in Fig. |7J The fraction of molecular to atomic gas in 
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Table 5. Ancillary data used in this paper. 



Wa vplpn trth 


Rpt;n1iiti on 

J VJ 1 LI L I VJ 1 1 




160 pm 


40" 


Spitzer MIPS a 


70 pm 


18" 


Spitzer MIPS a 


24 pm 


6" 


Spitzer MIPS a 


8 pm 


2" 


Spfeer IRAC b 


6570 A (Ha) 


2" (pixel size) 


KPNO c 


1500 A 


5" 


GALEX FUV d 



Table 6. Exponential scale length L in kpc for Spitzer MIPS 
and CO emission computed over the area mapped in CO. 





[0.5-3.5] kpc 


[3.5-7] kpc 


[0.5-7] kpc 


24 pm 


1.36 + 0.05 


1.02 + 0.07 


1.40 + 0.03 


70 pm 


1.51+0.07 


1.05 + 0.06 


1.48 + 0.04 


160 pm 


2.26 + 0.12 


1.30 + 0.06 


1.83 + 0.07 


CO 


2.0 + 0.2 


1.6 + 0.3 


1.9 + 0.1 



" iHinz et alj d2004l);lTabatabaei et al.l l l2007bh 

h IVerlevetal.l([2007h 

c iHoopes et al.l 12001): Gre enawa 

d iThilker et alj (l2005h . retrieved from GR5 public release of the 
MAST archive 
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Fig. 8. Top: Radial distribution of the Ho and Hi mass surface 
density. Bottom: Radial distribution of the molecular to atomic 
mass fraction. 



mass is shown in the bottom panel of Fig. [8] Given the constant 
atomic gas surface density, the fraction decreases following the 
molecular gas density from 1.2 to 0.015 from the center to the 
outer parts of the stellar disk. 

3.3.3. Star Formation Efficiency 

The star formation efficiency is usually defined as the ratio of the 
star formation rate over the mass of molecular gas available to 
form stars. 

SFE = SFR/M H2 

Defined in this way, the SFE is the inverse of the molecular 
gas depletion time. Estima t es from a variety of observations 
dBlitz & Rosolowskvl 120061: iKennicutj 119981: iHippelein et all 



120031) give a value of the SFR which varies from 0.3 to 
0.7 M ^yr -1 for the entire disk. A recent study by IVerlev et al.l 
(120091) using extinction corrected Ho-, FUV, and infrared tracers 
gives a value of 0.45 + 0.1 M0yr~' for total star formation rate. 
This is our adopted value. 

Our CO map does not cover the whole disk but covers re- 
spectively 63%, 61% and 59% of the IR emission at 24 /mi, 
70 pm and 160 pm (see bottom panel of Fig. [7]) so the SFR is 
~ 0.27 Mq yr -1 within the area covered by our CO map. With a 
Af(H 2 )/Ico(i-0) value of4xl0 20 cm~ 2 /(K km s~') and an average 
2-1/1-0 line ratio of 0.73 (see Sect. 13. 1 . ll for details), the molec- 
ular gas mass within the CO map is 1 .9x 10 8 Mq, With these val- 
ues, the SFE is 1 .4 x 10~ 9 yr~' or a molecular gas depletion time 
7 x 10 8 yr. An SFE of 1.6 x 10 9 yr _1 is a facto r ~ 3 higher than 
that found in large spiral galaxies (s ample in. lKennicuttlll998l: 
iMurgia et al.ll2002l: iLerov et a l. 2008) Earlier studies used lower 
Af(H2)/Ico(i-0) values and thus obtained even higher SFEs (or 
shorter H2 depletion times) as is shown in Table [4] 

In order to bring the SFE of M33 to the same level as in large 
spirals, the A^(H2)/Ico(i-0) factor would have to be about 3 times 
the value we use whereas previous studies used even smaller val- 
ues (I Wilson & Scovilldll990l:ICorbellill2003l:lHeyeret al.ll2004t 
lEngargiola et al.l 12003b Rosolowsky et al r i2007b Garda n et al.l 
|2007|) . Through Virial mass estimates and 13 ^ 12 CO line ratios, 
Braine et al (submitted to A&A) support a value of about 
W(H 2 )/I C o(i-o) =4xl0 20 cm^/CKkms-^inthediskofMSS, 
possibly closer to the Galactic value in the inner disk. 

Figure |9]presents the radial distributions of the SFE as cal- 
culated from the Ha, 24, and 70 pm emission. The units are 
arbitrary as we are interested in showing the radial trend. The 
SFR were taken as being directly proportional to the Ha and IR 
luminosities. Ha and FIR emission tend to show opposite biases 
- the Ha suffers extinction where the FIR is strong and the FIR 
is weak where little dust is present but the Ha is unaffected. If 
one were to compute the star formation efficiency taking into ac- 
count the total gas mass, the SFE would drop dramatically with 
increasing radius as opposed to the SFE computed using only 
molecular gas that does not appear to vary with radius. Table |4] 
summarizes the masses, star formation rates and depletion times 
for this study and previous surveys of M33. 



3.4. The ISM and star formation tracers 

In this section, we compare the CO emission with tracers of star 
formation. We present a series of figures ([T0lto[T4l) comparing 
the CO emission, as contours, overlaid on images of the 8, 24 
and 70 pm FIR emission as well as the Ha and FUV emission 
(Blow ups of these 5 figures are available in the online version of 
this article). The 160 pm image is at lower resolution and brings 
little morphological information not present in the 70 pm image. 
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Fig. 9. Radial distribution of the Star Formation Efficiency in 
arbitrary units, all of the curves are normalized to the central 
value and shifted by 1 decade along the vertical axis. 



Table 7. Comparison with previous surveys of CO emission 
exponential scale lengths L in kpc. 



Study 


Range 


Scalelength 




(kpc) 


(kpc) 


This paper 


0-7 


1.9 + 0.1 


Gardan et al. (7007) 


2-6 


1.4 + 0.1 


Engargiola et al. (2003) 


0-7 


1.4 + 0.1 


Corbelli (2003) 


0-6 


2.5 



While the dynamic range in the CO image is much lower 
than in the images of star formation tracers, the correspondence 
is globally excellent. Figure [16] shows the scatter plot between 
the CO integrated intensity and 70 /mi emission. 

Figures [TOl4T4"l show that in almost all cases the CO emission 
is found where star formation is detected and in general the CO 
emission is stronger as the level of SF increases. The weaker 
CO contours outline the regions with a remarkable precision (see 
online figures 1 to 12 for the blowups). A partial exception is 
found near (a = l h 34"W4, 6 = +30°49'06", J2000) where 
the CO emission is extremely strong while the SF tracers show 
rather weak emission. The CO spectrum in the direction of that 
particular cloud is shown in Fig. [17] 

A true exception to the FIR-HI-CO correlation is the inter- 
arm cloud at (a = l h 34 m 16 s :7, 5 = +30°59'3", J2000) dubbed 
"Lonely Cloud" bv lGardan et alJ (120071) . It is found in a region 
with very little if any star formation and relatively weak H i emis- 
sion, despite fairly strong CO (and even 13 CO) emission. When 
H i brightness thresholds higher than 10 K are used to create the 
masked zero-moment CO maps, the Lonely Cloud disappears. 
While this cloud was discovered in the first set of observations, 
it remains a unique object as the observations of the central and 
southern regions, de spite considerably lower noise levels than in 
iGardan et al.l (|2007), have not revealed similar molecular clouds. 



WBPW 




Fig. 10. Galex FUV image with Ico main beam contours of 
1 (white), 2, 4 (black), 8 (blue) Kkms -1 . The beam size is 
shown as a white dot in the lower left corner. Blow ups are avail- 
able in the online version. 



In the outer parts, the FIR emission is typically quite incon- 
spicuous, even where CO has been detected, but the Ha emis- 
sion is often quite strong. In the central regions, the alignment of 
the emission at the different wavelengths is nearly perfect, espe- 
cially between CO and FIR. Further out, however, notable offsets 
become more common although there is a clear association be- 
tween the CO and FIR/Ha emitting regions. The offsets are less 
than 100 pc, particularly in comparison with the FIR. 

The CO emission follows Hi-defined arm-like structures 
very closely, completely avoiding the Hi holes. As can be seen 
in Fig. [15] and corresponding blow ups in the online version of 
the paper, the overall correlation between H i and CO emission is 
very good. CO is generally found in Hi-rich regions although not 
all of the H i peaks show CO emission and not all CO emission 
falls on H i. It is not clear what the source of the 70 /mi emission 
is in for example the Hi (and CO) hole around (a = l h 34 m 25, 
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Fig. 11. KPNO Ha image with Ico main beam contours of Fig. 12. Spitzer 8 /mi image with Ico main beam contours of 
1 (white), 2, 4 (black), 8 (blue) Kkms -1 . The beam size is 1 (white), 2, 4 (black), 8 (blue) Kkrns -1 . The beam size is 
shown as a white dot in the lower left corner. Blow ups are avail- shown as a white dot in the lower left corner. Blow ups are avail- 
able in the online version. able in the online version. 



8 = +30°35'00", J2000) as the amount of neutral gas is very 
low. The emission in this region is lower in the 24 and 8 /urn 
wavebands, so perhaps the 70 micron emission appears stronger 
than in reality due to the poorer spatial resolution. The upcoming 
Herschel HERM33ES observations should answer this. 

3.5. Examples of interesting spectra 

The spectra shown in Figures [TTJ [18] and [19] are just a few 
among a huge number but are intended to show "special" po- 
sitions. Fig [17] shows a very strong CO profile towards a re- 
gion without strong st ar formation (note good agreement with 
lEngargiola et al. 20031). The line is very strong with a line width 
of about 11 kms~' at half-intensity. The other spectra shown 
are interesting because the CO lines are detected very far from 
the center of M33, despite the low metallicity and surface mass 



density. Fig. [18] shows the second molecular cloud detected at 
or beyond R25 in M33. A secondary peak offset in velocity by 
about 8 km s _1 may be present. Fig[l9]shows further spectra near 
R25 which clearly is composed of emission from two clouds, the 
cloud to the SE (of the figure) at -260 km s _1 and the one to the 
NW at about -273 km s _1 . The spectral resolution of the data is 
barely sufficient to measure the linewidths of the clouds, which 
are typically less than 5 km s _1 at half p ower. This was also t rue 
of the spectra of distant clouds shown in lBraine et al 
they are much narrower than in the H2-rich inner disk. It is not 
clear whether the broad linewidths in the inner disk are because 
the GMCs themselves are much more massive or due to a super- 
position of clouds (or a single GMC in a region with significant 
diffuse emission). 

The spectra shown also illustrate the dynamic range of the 
data acquired, roughly 100 in integrated intensity. 
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Fig. 13. Spitzer 24 //m image with Ico main beam contours of 
1 (white), 2, 4 (black), 8 (blue) Kkms~'. The beam size is 
shown as a white dot in the lower left corner. Blow ups are avail- 
able in the online version. 
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Fig. 14. Spitzer 70 fim image with Ico main beam contours of 
1 (white), 2, 4 (black), 8 (blue) Kkms -1 . The beam size is 
shown as a white dot in the lower left corner. Blow ups are avail- 
able in the online version. 



4. Conclusions 

This work presents high-resolution maps of the atomic and 
molecular gas in the disk of M33 via observations of the 
Hi 21 cm and CO(2-l) lines. The whole disk out to 8.5 kpc 
is covered in the H i line and about 60% of the emission in CO. 

Assuming the iV(H2)/Ico(t-0) factor to be twice that of the 
Galaxy, because of the sub-solar metallicity of M33, we esti- 
mate a molecular gas mass of 3.3 x 10 8 Mq roughly 20% of the 
1.4 x 10 9 M0 detected in the inner 8.5 kpc in Hi. Azimuthally 
averaging, the H i surface density is close to constant with radius 
but the H2 decreases exponentially with a scale length of 1.9kpc. 
The H2/H1 mass ratio decreases from about unity to 1%. 

The correspondence between the peaks and holes in the dis- 
tributions of molecular and atomic gas is excellent and follows 
the peaks and troughs in the FIR, MIR, and Ha images. The SFE 



is approximately constant with radius, suggesting that molecu- 
lar gas is transformed into stars at a similar rate (assuming a 
similar IMF) at all galactocentric radii. However, the SFE in 
the small, gas-rich, low-metallicity, blue spiral M33 appears 2-4 
times higher than what is observed in large nearby spirals. 

The sensitivity of the survey is such that CO emission is de- 
tected far out in the disk of M33 although few clouds are present 
and the lines are much weaker in intensity. 

Acknowledgements. We thank the IRAM staff in Granada for their help with the 
observations. 
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Fig. 17. Spectrum of a bright cloud at position (a = l h 34 m 09"4, 
Fig. 15. VLA H i zeroth moment image with I co main beam con- $ = +30 °49'06", J2000), unlike other strong peak of CO inte- 
tours of 1 (green), 2 (red), 4(black), 8 (blue) K km s .The grated intensityi it lacks strong emission in star formation tracers 
beam size is shown as a white dot in the lower left corner. Note ^ see Sectf34| 
how the CO emission follows the bright H i and outlines the H i 
holes. Blow ups are available in the online version. 
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Fig. 18. Spectra of the outer disk cloud (a = l h 34 m 50:8, 6 = 
+31°08'26", J2000) shown with 6" spacing between panels. If 
the secondary peak at -266 km s _1 , seen in the second spectrum 
from the bottom in the middle column, is real then the cloud is 
small as it is only seen within a single beam. 




0.05 


g 

% 0.05 

I ' 



0.05 - 



.1 1 1 1 1 1 1 1 1 1 1 1 1_ 

j 1 1 1 1 1 1 1 1 1 1 1 \ 


_i 1 1 1 1 1 1 1 1 1 1 1 1_ 
) 1 1 1 1 1 1 1 1 1 | f 


) 1 1 1 1 


i 1 1 i i 1 1 i 
I I | I I I | I 


.i i i i i i 1 1 i i i i i, 
j I I | I I I | I I I | f 












"i 1 1 yi i 1 1 1 1 1 1 1" 
j 1 1 1 1 1 1 1 1 1 1 1 i 


"i i i 1 i i i 1 i i i 1 r 
j 1 1 1 1 1 1 1 1 1 1 1 r 


■ 1 1 1 1 1 1 1 1 1 1 1 1 

j 1 1 1 1 1 1 1 1 1 1 1 1 


"i i i | i i i | i i i | r 

j 1 1 1 1 1 1 1 1 1 1 1 r 


■ , , i , , , i , , , i ■■ 


-, , , i , , , i , , , i ,- 


; , , I , , , i , , , i , 


-. , , i . , . i , , , i ,- 



-200-360-240 - EBO-E60-210 -2B0-EEO-E4D -E8O-260-E4O 
Velocity (km/s) 

Fig. 19. Spectra of two more outer disk clouds at different veloc- 
ities. The lower left panel is at (or = l h 34 m 21!7,5 = +31°04'13", 
J2000) and the panel spacing is 6". 
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Appendix A: Online Figures 




i^iW Htf 1 GO 1 

Fig. A.l. Northern part of KPNO Ha image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s . The beam 
size is shown as a white dot in the lower left corner. 



P. Gratier et al.: Molecular and Atomic Gas in the Local Group Galaxy M33 17 




I^jfl" BOP 



Fig. A.2. Center part of KPNO Ha image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam size 
is shown as a white dot in the lower left corner. 
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Fig. A.3. Southern part of KPNO Ha image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 
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Fig. A.4. Northern part of Spitzer 8 fim image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 
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Fig.A.5. Center part of Spitzer 8 fim image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) Kkm s _1 . The beam 
size is shown as a white dot in the lower left corner. 
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Fig. A.6. Southern part of Spitzer 8 fim image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 




Fig. A.7. Northern part of Spitzer 24 fim image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s . The beam 
size is shown as a white dot in the lower left corner. 
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Fig. A.8. Center part of Spitzer 24 fim image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 
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Fig.A.10. Northern part of Spitzer 70 /im image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) Kkm s The 
beam size is shown as a white dot in the lower left corner. 
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Fig. A.ll. Center part of Spitzer 70 yum image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 
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Fig. A.12. Southern part of Spitzer 70 /im image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) Kkm s The 
beam size is shown as a white dot in the lower left corner. 




Fig. A.13. Northern part of VLA H i 2 1 cm integrated intensity image with Ico main beam contours of 1 (green), 2 (red), 4(black), 
8 (blue) K km s _1 . The beam size is shown as a white dot in the lower left corner. 
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Fig. A.14. Center part of VLA Hi 21 cm integrated intensity image with Ico main beam contours of 1 (green), 2 (red), 4(black), 
8 (blue) K km s -1 . The beam size is shown as a white dot in the lower left corner. 
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Fig. A.15. Southern part of VLA H i 21 cm integrated intensity image with Ico main beam contours of 1 (green), 2 (red), 4(black), 
8 (blue) K km s _1 . The beam size is shown as a white dot in the lower left corner. 
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Fig. 9. Southern part of Spitzer 24 yum image with Ico main beam contours of 1 (white), 2, 4 (black), 8 (blue) K km s 1 . The beam 
size is shown as a white dot in the lower left corner. 



